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SEEING DOUBLE: STRONG GRAVITATIONAL LENSING OF HIGH-REDSHIFT SUPERNOVAE 
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ABSTRACT 

With the advent of large, deep surveys, the observation of a strongly gravitationally lensed supernova 
becomes increasingly likely. High-redshift surveys continue apace, with a handful of type la supernovae 
observed to date at redshifts of one or greater. In addition, a satellite (the Supernova/ Acceleration 
Pro6e[SNAP]) has been proposed dedicated to observing thousands of supernovae per year out to a 
redshift of 1.7. Although it is exceedingly unlikely that we will see a multiply- imaged supernova from 
ongoing surveys, we find that SNAP would observe at least eight such events per year. Since having 
a standard candle is inessential to most lensing studies, SNAP's large sample of type II supernovae 
contributes to this rate. Each case of strong lensing allows for a precise determination of time delays, 
image separations, and relative image magnifications, and the SNAP strong-lensing database will offer 
measures of ^^A, and Hq, independent of SNAP's primary goal of establishing the distance-redshift 
relation. These systems also constrain models for the matter density profiles of galaxies and clusters. 
Furthermore, lensed type la supernovae afford the opportunity to break the mass-sheet degeneracy found 
in many lensing measurements. 

Subject headings: cosmology: observations — cosmology: theory — gravitational lensing — supernovae: 
general 



1. INTRODUCTION 

Type la supernovae (SNe) have become an extremely 
powerful cosmological probe. They appear to be very 
good (cahbrated) standard candles (Riess, Press, & Kir- 
shner 1996a; Branch 1998), and because of this they en- 
able a direct measurement of the distance-redshift rela- 
tion, and hence the expansion history of the universe (Col- 
gate 1979; Goobar & Perlmutter 1995). At present on the 
order of 100 high-redshift SNe have been observed, and 
surveys are adding dozens of type la's a year at increas- 
ingly high redshifts. Results from these surveys indicate 
that the universe's expansion is accelerating (Riess et al. 
1998; Perlmutter et al. 1999), and this discovery has led 
to a tremendous amount of excitement in the cosmological 
community. In particular, there is great interest in the Su- 
perNova/ Acceleration Probe (SNAP), a satellite dedicated 
to searching for and observing high-z SNe.^ The primary 
objective of the mission would be to measure cosmolog- 
ical parameters to unprecedented accuracy, enabling, for 
example, the distinction between a cosmological constant 
and a "quintessence" component (Huterer & Turner 1999). 
The satellite would achieve this by observing over 2000 
type la SNe per year, getting a complete light curve, mul- 
tiple colors, and a spectrum at maximum light for each 
one, and thereby measuring the distance-redshift relation 
out to z = 1.7 to unprecedented precision. 

Gravitational lensing can influence the observed bright- 
ness of high-z SN data. The impact of small shifts in the 
brightnesses of SNe due to weak lensing has been exten- 
sively explored, both as a potential source of noise (Kan- 
towski, Vaughan, & Branch 1995; Frieman 1997; Wamb- 
sganss et al. 1997; Holz 1998), and as a source of sci- 
ence (Metcalf & Silk 1999; Holz 1999; Seljak & Holz 1999; 
Metcalf 1999). In this paper we focus on the possibil- 

^ see http://snap.lbl.gov for more information. 



ity for strong lensing of SNe, leading to multiple images 
of a given source. While this possibility has also been 
discussed (Schneider & Wagoner 1987; Lindcr, Schneider, 
& Wagoner 1988; Kovner & Paczyhski 1988; Porciani & 
Madau 2000), we consider strong lensing in the age of 
large-scale, high-z SN surveys (see also Wang (2000)). Al- 
though current surveys are unlikely to produce a case of 
multiple imaging, we find that SNAP would be expected 
to see numerous such events. In addition, SNAP would be 
well suited to do extensive characterization of the events, 
bringing the very real possibility of science from multiply- 
imaged SNe in the foreseeable future. 

2. LIKELIHOOD OF MULTIPLY-IMAGED SUPERNOVAE 

We commence by estimating the likelihood of strong 
lensing. An empirical value comes from the rate of ob- 
served strong-lensing in optical quasar or radio source 
surveys. In an attempt to control selection effects, the 
Jodrell-Bank VLA Astrometric Survey (JVAS) and the 
Cosmic Lens All-Sky Survey (CLASS) have conducted a 
statistically complete survey of over 15,000 flat-spectrum 
radio sources (Brown & Myers 2000; Myers et al. 2001). 
With analysis almost complete, they have found 18 lensed 
sources, with 14 having source redshifts and image time de- 
lays that fall within the SNAP detectability range. Taken 
at face value (e.g. neglecting selection effects and popu- 
lation evolution), this would indicate that approximately 
0.1% of SNe seen by SNAP would be expected to be 
strongly lensed, with image separations on the order of 
one arcsec and time delays between images ranging from 
a few days to over a year. As this statistic is based on a 
small number of strong lensing events, it is also instructive 
to look at theoretical calculations of the lensing rate. With 
this in mind we utilize the "stochastic universe method" 
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(SUM) to determine the probability for multiple imag- 
ing (Holz & Wald 1998). A key feature of this approach 
is that it accommodates lensing from all matter along 
the line of sight — possibly including the effects from thou- 
sands of galaxies and in addition requires no assumptions 
regarding the luminosity-distance relation (Holz & Wald 
1998; Bergstrom et al. 2000). It is to be emphasized, how- 
ever, that the strong lensing results do not depend sensi- 
tively on the details of the methods, with differing analytic 
and numerical results being in good agreement (Turner, 
Ostriker, & Gott 1984; Holz, Miller, & Quashnock 1999; 
Metcalf & Silk 1999). 

In Fig. 1 we show the probability of multiple-imaging 
of a given source (also referred to as the optical depth 
to strong lensing, r) as a function of the source red- 
shift, for a range of cosmologies. We take a conserva- 
tive value of the dimensionless lensing efficacy parameter, 
F = 0.025 (Turner ct al. 1984). The results arc relatively 
insensitive to details of the mass function, and calcula- 
tions with a PS distribution (Press & Schechter 1974) of 
galaxy masses yields similar results. Although taking all 
of the matter in the Universe to be in isothermal spheres 
is a great oversimplification, it is to be emphasized that 
effective modeling of many strong-lensing systems is ac- 
complished with isothermal mass distributions (Koopmans 
& Fassnacht 1999; Witt, Mao, & Keeton 2000; Cohn et 
al. 2001).2 Although NFW profiles (Navarro, Frenk, & 
White 1996) arc expected to be a better approximation 
to the mass distribution of massive halos, baryons are ex- 
pected to isothermalize the matter distributions for halos 
of galaxy mass and below (Kochanek & White 2001). Por- 
ciani & Madau (2000) have undertaken a study of lensing 
in cold dark matter universes with a PS distribution of 
galaxy masses, where NFW profiles are utilized for mas- 
sive halos (m > 3.5 x 10^^ M0), and singular isothermal 
spherical profiles are used otherwise. Our extremely sim- 
plified model agrees with their results to within ~ 20%. 
Our results are also consistent with the CLASS data. 

Data in cosmology in recent years has been pointing 
to a "cosmic concordance" model (Bahcall et al. 1999), 
with n„, ^ 1/3, = 2/3, and Hq = 70 kmsec^^ Mpc"\ 
and we restrict our attention to this model in what fol- 
lows. The likelihood results in Fig. 1 imply that 0.05% of 
sources at rcdshift z = 1 would be expected to be multiply 
imaged on arcsecond scales. At present roughly 100 high- 2: 
SNe have been observed (mostly at 2: « 0.5), and thus it is 
quite unlikely that any of them are multiply imaged. The 
situation is more bleak than the numbers indicate, how- 
ever, as even if one of the observed high-z SNe happens 
to be multiply-imaged, it is exceedingly unlikely that we 
would stumble upon successive images. 

The situation is much brighter if SNAP flies. SNAP's 
"optical" (0.35-1.0 ^m) imager would observe 3,800 type 
la SNe per year, with detailed follow-up (restframc B- 
band photometry and spectra) of a subsample of 2,400 
(100 at z > 1.2, and all at z < 1.2). By checking 
whether prior images have appeared nearby (within 15") 
when selecting high-z SNe for follow-up, SNAP can assure 
that every multiply-imaged type la SN has at least one 
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Fig. 1. — The optical depth, t, to multiple imaging, as a function 
of redshift, z. Results for a range of cosmological parameters are 
shown, for a lensing efficacy of F = 0.025. 

image with a redshift and absolute magnification deter- 
mination. The ratio of the optical fluxes then allows a 
determination of the absolute magnification of all other 
images. Convolving the optical depth curve of Fig. 1 with 
the expected redshift distribution of SNAP SNe (relatively 
flat, with 1, 500 type la's/year at 1.2 < z < 1.7) yields an 
expectation of 2 multiply imaged la SNe per year. 

SNAP would provide a homogeneously selected 28.5 
(restframe U or B-band) magnitude-limited sample of all 
SNe which occur during its observation period in its 20 
flxed 1 deg^ fields of view, with two fields being addition- 
ally imaged down to a limiting magnitude of 30. Since the 
gap between deep images of each field will always be less 
than eight days, no SNe will fall through the observing 
cracks. AlthoTigh the SNAP design is still in flux, we will 
take as one of its core reqtiirements that all type la SNe 
with 2 < 1.7 will be discovered at least 3.8 magnitudes 
below peak brightness (and thus within ~ 2 rest frame 
days of the explosion). Therefore lensed images which 
have > 3% of the (unlensed) source flux will be observed. 
From our simplified SUM calculations, we find that 94% 
of strongly-lensed events are at least this bright. SNAP 
will see every type la SN out to redshift 1.7 in its deep 
fields, and will observe the vast majority of strong lensing 
events of these sources, thereby minimizing magnification 
bias in the strong lensing sample. 

SNAP will also see a very large sample of type II SNe 
(though only a tiny fraction of these will have detailed 
follow-ups). While these are significantly dimmer (peak 
value 2.3 magnitudes below the la peak (in restframe B- 
band), with a dispersion of 1.3 mag (Gilliland, Nugent, & 
Phillips 1999)), they also occur at much higher rates (3 
times more frequent at present day. rising to 5-10 times 
more frequent at redshift 1 2 (Sullivan et al. 2000)). Tak- 
ing a conservative value of 5 times the la rate, we find 
that SNAP would have 10 multiply-imaged type II SNe in 
its sample per year. SNAP will catch the average SN II 
event only 1.5 magnitudes below peak, and so will miss 
images demagnified to less than 25% of the unlensed flux. 
Convolving the large intrinsic dispersion in type II peak 
brightness with the distribution of image pair magnifica- 
tions (from our SUM calculations), we find that 60% of 
multiply-imaged type lis will be observable, yielding an 



^ Singular isothermal spheres always produce double images, while roughly half of observed strong lensing systems are quads (Rusin & Tegmark 
2001). Adding background shear, or considering singular isothermal ellipsoid lenses, produces the requisite quad images. Such generalizations 
can be expected to raise the optical depth curves in Fig. 1 (Schneider, Ehlers, & Falco 1992). 
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(^10%) fraction of halo matter is in macroscopic form (e.g. 
MACHOs, black holes), the lensing probability increases 
dramatically due to frequent microlensing (Schneider & 
Wagoner 1987; Ranch 1991). At cosmological distances, 
lensing by compact objects with masses below 10^ Mq 
produces image separations of less than a few milliarc- 
seconds and time delays of less than a minute, indicat- 
ing that multiple imaging would no longer be directly dis- 
cernible. The overall amplification of the combined im- 
age would still be noticeable down to lens mass scales of 
~ 10""' Mq, below which the SN can no longer be treated 
as a point source (linear extent ~ 10^^ cm corresponds to 
an angular size on the order of the Einstein angle of the 
lens), and lensing ceases to have an effect. In addition, 
if there is relative motion between the lens and the line 
of sight to the source, the amplification becomes time de- 
pendent. The lensing timescale at cosmological distances 
goes as 5t - 18 days (10^ cmsec-7^^)y(A'-^/lO"'^M0), 
with V the expansion velocity of the SN, and m the mass 
of the compact lens. Therefore lenses in the mass range 
10"'' Mq <m < 10-3 Mq can produce a MACHO-type 
signal superposed on the SN light curve. Microlensing by 
compact objects with m > 10^^ Mq will manifest itself as 
an overall amplification shift of the (standard candle) light 
curve, inducing further scatter in the SN Hubble diagram 
(Minty, Heavens, & Hawkins 2001). 

4. INDEPENDENT MEASURE OF COSMOLOGICAL 

CONSTITUTION 

Over the course of a few years, SNAP would provide a 
uniformly selected survey of many thousands of SNe, offer- 
ing the possibility of a variety of statistical tests of cosmol- 
ogy. For example, the different curves in Fig. 1 should be 
distinguishable observationally — three years' data would 
difi'erentiate between 11a = 2/3 and I^a = at better than 
3iT. That the SNe are standard candles is completely ir- 
relevant to determining cosmology via lensing likelihood, 
and thus strong lensing provides a powerful independent 
probe, and consistency check, for the SN measurement 
of the distance-redshift relation (the primary mission of 
SNAP). 

In addition, the strong lensing likelihood depends upon 
the dumpiness of the dark matter, and thereby allows 
the distinction between microscopic and macroscopic dark 
matter (Linder et al. 1988). Furthermore, Wyithe et al. 
(2001) have suggested that the optical depth to multiple 
imaging could be used to distinguish between interacting 
and non-interacting dark matter. A complementary cos- 
mological measurement is offered by the distribution of im- 
age separations and lensing image morphologies (Keeton, 
Kochanek, & Seljak 1997; Rusin & Ma 2001; Kochanek & 
White 2001). These approaches are less sensitive to the 
cosmological parameters, but quite sensitive to the mass 
distribution within halos. 

5. INDEPENDENT MEASURE OF COSMOLOGICAL SCALE 

Multiply-imaged systems allow for the determination of 
the Hubble constant, Hq, through the measurement of 

^ Extinction along different lines of sight could yield a chromatic difference between multiple images of a given SN. However, because of the 
extensive phenomenological knowledge of type la SN light curves and colors (Riess, Press, & Kirshner 1996a,b), the extinction in a given image 
is quantifiable to within < 0.04 mag, sharply reducing the likelihood of false dissociations. Having multiple images sampling different lines of 
sight through a galaxy or cluster could provide a rudimentary map of the distribution of dust in the lens. 



expected total of 6 events per year. For these systems, 
snap's optical imager light curve will provide a precise 
determination of the image time delays, separations, and 
relative fluxes. As one does not require the standard can- 
dle luminosities of Icnsed SNe to do strong-lcnsing science 
(see Sec. 4 and 5), strongly lensed type II systems are a 
valuable contribution to the lens database. 

Combining the rate for SNe la and II, we conservatively 
estimate that SNAP would detect eight multiply-imaged 
SNe per year. Lot us therefore assume that SNAP will fly, 
and explore the science that can be done with multiply- 
imaged SN systems. 

3. DETECTION CRITERIA 

A multiply- imaged SN would appear as successive SNe, 
closely clustered on the sky (< 15"). As gravitational lens- 
ing is achromatic, the light curves of the "different" SNe 
would be identical in all bands, modulo an overall amplifi- 
cation shift between the different images. ^ The probability 
of having a second SN appear in a given galaxy within one 
year is roughly equivalent to the multiple-imaging rate, 
and therefore well-sampled photometry of the light curves 
(to enable comparison of rise/fall parameters), or searches 
for lenses along the line of sight, will be required to posi- 
tively identify lensing cases. 

Characteristics of the SNAP mission impose a number 
of selection constraints on observable multiply-imaged SNe 
systems. The distribution of time delays between SN im- 
ages can be expected to mimic those in multiply-imaged 
radio source and quasar surveys. Prom the theoretical dis- 
tributions of time delays obtained via the SUM calcula- 
tions, we find that 78% of strong lensing events have time 
delays such that they'd be seen in one year's observation, 
rising to 93% with three years of observation. 

The angular resolution of SNAP is to be O'.'l, which 
allows it to directly resolve the vast majority of image 
separations. An important possible systematic error for 
ground-based optical surveys of multiply-imaged systems 
is that images appearing near the center of the lens can be 
obscured by the lens (Porciani & Madau 2000). Because 
of snap's excellent angular resolution, and because base- 
line images of the galaxy without the supernovae will be 
available (and hence can be subtracted off), a SN image 
would be identifiable all the way into the core of a galaxy. 

Even if the angular separation of SN images is well be- 
low the resolution limit of the SNAP telescope (which is 
exceedingly unlikely for macrolensing cases), because the 
SNe are transitory we may nonetheless be able to see and 
distinguish images. As long as the light curves do not 
significantly overlap in time, proximity (and blending) of 
images is no impediment to observation. In addition, it is 
straightforward to associate even very widely spaced im- 
ages, as we expect to see all SNe within SNAP's large 
(1 deg^) field of view (limited to systems with time delays 
which fall within the SNAP observational period). 

The likelihood results are quite sensitive to the form 
of the dark matter. Thus far wc have been considering 
smoothly distributed isothermal halos filled with micro- 
scopic dark matter (e.g. axions, WIMPS). If a significant 
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time delays between images (Refsdal 1964). This method 
has been developed for over a decade in multiply-imaged 
quasar systems, with varying success (Koopmans & Fass- 
nacht 1999). One of the limitations of such systems is 
that the time delays can often be quite difficult to deter- 
mine (Kundic et al. 1997). In the case of multiply-imaged 
SNe the determination of the time delay will be completely 
straightforward, as the time of peak luminosity of a given 
SN will be determined to an accuracy of better than a 
day as an integral part of the SNAP observing program. 
It is also often difficult to determine the flux ratio of the 
images with good accuracy (due to time variability of the 
source, extinction, etc.; see Mao & Schneider (1998)). As 
SNe light curves have been extensively characterized, the 
flux ratios will be well determined. 

It is often the modeling of the lens that limits the preci- 
sion of the determination of Hq (Keeton et al. 2000; Cohn 
ct al. 2001). A major impediment to constraining lens 
profiles is the limited number of constraints a lens system 
provides (at best, a handful of image separations, rela- 
tive brightnesses, and time delays). Attempts to increase 
the number of constraints generally rely on utilizing lensed 
images of the extended background galaxy (arcs, Einstein 
rings, etc.; see Keeton et al. (2000); Kochanek, Keeton, & 
McLeod (2001)). For quasars these attempts can be hin- 
dered by the large brightness contrast between the quasar 
source and the rest of the host galaxy. In the case of SNe, 
however, the images conveniently turn themselves off, al- 
lowing an uncontaminated view of the host. To establish 
a baseline for the light "contamination" to the SN curve 
from the host galaxies, SNAP will make deep images of 
the hosts both before the SNe appear and after they fade 
away. Subsequent very deep observations with a variety 
of different instruments (e.g. NGST, CELT+AO) at a 
range of wavelengths could also be made. These images of 
the lensed host galaxy would provide important additional 
constraints to the lens model. It is again to be emphasized 
that the measure of Hq with time delays is independent 
of the standard-candleness of the SNe, and is thus com- 
pletely independent of Cepheid-based measurements of the 
Hubble constant with type la SNe (Jha et al. 1999). 



6. BREAKING THE MASS SHEET DEGENERACY 

Although the science in the preceding sections does not 
depend on the SNe being standard candles, this property 
can have important consequences for lensing studies. A 
drawback of many lensing analyses done to date is that 
they suffer from a mass-sheet degeneracy: a uniform sheet 
of matter anywhere between the source and observer will 
generally remain undetected. Because type la (and to a 
lesser degree type II) SNe are standard candles, the ab- 
solute magnification can be determined, yielding an ad- 
ditional constraint (by changing the relative fluxes of all 
images into absolute fluxes), and breaking this degener- 
acy (Kolatt & Bartelmann 1998). 

An additional feature of SNAP is that, in the course 
of returning to the same flelds week after week looking 
for new SNe, it will make very high quality weak lens- 
ing maps. Coupling information from both weak and 
strong lensing studies of a given fleld may yield additional 
insights — for example, a mass concentration which causes 
a large-separation multiple-imaging event may also be de- 
tectable in the weak lensing map of the same field. Any 
SNe (strongly lensed or not) superposed on a weak lensing 
map will (in principle) allow for a breaking of the mass 
sheet degeneracy. 

7. CONCLUSIONS 

The observation of multiply- imaged SNe would afford 
independent measurements of cosmological parameters. 
Although current high-^ SN surveys are unlikely to ob- 
serve a multiply-imaged SN system, SNAP would be ex- 
pected to sec at least eight such events per year. These 
systems would be a boon to science, independent of and 
in addition to the primary goals of SNAP. 
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NSF under grant PHY99-07949 to the ITP. 



REFERENCES 



Bahcall, N., Ostriker, J. R, Perlmutter, S., & Steinhardt, P. 1999, 

Science, 284, 1481 
Bergstrom, L., Goliath, M., Goobar, A., & Mortsell, E. 2000, A&A, 

358, 13 

Branch, D. 1998, ARA&A, 36, 17 

Browne, I. W. A., & Myers, S. T. 2000, lAU Circ, 201, 47 
Cohn, J. D., Kochanek, C. S., McLeod, B. A., & Keeton, C. R. 

2001, ApJ, 554, 1216 
Colgate, S. 1979, ApJ, 232, 404 

Frieman, J. A. 1997, Comments Astrophys., 18, 323 

Gilliland, R. L., Nugent, P. E., & Phillips, M. M. 1999, ApJ, 521, 30 

Goobar, A. & Perlmutter, S. 1995, ApJ, 450, 14 

Holz, D. E. 1998, ApJ, 506, LI 

Holz, D. E. 1999, Nature, 400, 819 

Holz, D. E., Miller, M. C, & Quashnock, J. M. 1999, ApJ, 510, 54 
Holz, D. E., & Wald, R. M. 1998, Phys. Rev. D, 58, 063501 
Huterer, D., & Turner, M. S. 1999, Phys. Rev. D, 60, 081301 

Jha, S. et al. 1999, ApJS, 125, 73 

Kantowski, R., Vaughan, T., & Branch, D. 1995, ApJ, 447, 35 
Keeton, C. R., Kochanek, C. S., & Scljak, U. 1997, ApJ, 482, 604 
Keeton, C. R., et al. 2000, ApJ, 542, 74 

Kochanek, C. S., Keeton, C. R., & McLeod, B. A. 2001, ApJ, 547, 
50 



Kochanek, C. S. & White, M. 2001, astro- ph/0102334 
Kolatt, T. S., & Bartelmann, M. 1998, MNRAS, 296, 763 
Koopmans, L. V. E., & Fassnacht, C. D. 1999, ApJ, 527, 513 
Kovner, I., & Paczyhski, B. 1988, ApJ, 335, L9 
Kundic et al. 1997, ApJ, 482, 75 

Linder, E. V., Schneider, P., & Wagoner, R. V. 1988, ApJ, 324, 786 
Mao, S., & Schneider, P. 1998, MNRAS, 295, 587 
Metcalf, R. B. 1999, MNRAS, 305, 746 
Metcalf, R. B., & Silk, J. 1999, ApJ, 519, LI 

Minty, E.M., Heavens, A.F., & Hawkins, M.R.S. 2001, MNRAS, 

submitted (astro-ph/0104221) 
Myers, S.T. et al. 200i, in preparation; see the CLASS homepage 

http : //info . aoc .nrao . edu/'smyers/class .html 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 
Perlmutter, S. et al. 1999, ApJ, 517, 565 
Porciani, C, & Madau, P. 2000, ApJ, 532, 679 
Press, W. H., & Scheehter, P. 1974, ApJ, 187, 425 
Rauch, K. P. 1991. ApJ. 374, 83 
Refsdal, S. 1964, MNRAS, 128, 307 

Riess, A., Press, W. H., & Kirshner, R. P., 1996, ApJ, 473, 88 
Riess, A., Press, W. H., & Kirshner, R. P., 1996, ApJ, 473, 588 
Riess, A., et al. 1998, AJ, 116, 1009 
Rusin, D., & Ma, C.-P. 2001, ApJ, 549, L33 



Strong lensing of supernovae 



5 



Rusin, D., & Tegmark, M. 2001, ApJ, 553, 709 

Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses 

(Berlin: Springer) 
Schneider, P., & Wagoner, R. V. 1987, ApJ, 314, 154 
Seljak, U., & Holz, D. E. 1999, A&A, 351, LIO 
SuUivan, M., ElUs, R., Nugent, P., Smail, I., & Madau, P. 2000, 

MNRAS319, 549 



Turner, E.L., Ostriker, J. P., & Gott, J. R. 1984, ApJ, 284, 1 
Wambsganss, J., Gen, R., Xu, G., & Ostriker, J. P. 1997, ApJ, 475, 
L81 

Wang, Y. 2000, ApJ, 531, 676 

Witt, H. J., Mao, S., & Keeton, G. R. 2000, ApJ, 544, 98 
Wyithe, J. S. B., Turner, E. L., & Spergel, D. N. 2001, ApJ, 555, 
504 



